Fermiology of a ID heavy electron metal 
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We present a Quantum Monte Carlo (QMC) study of the ID Kondo lattice at non-integer filling, 
i.e. a one-dimensional version of a Heavy Fermion metal. For this special system the minus-sign 
problem turns out to be greatly reduced, so that accurate QMC simulations for temperatures as 
low as 1% of the conduction electron bandwidth are feasible. The single particle Green's function 
shows and intricate network of low-energy bands at low temperature, with a Fermi surface volume 
that comprises both c and /-electrons. As the temperature increases the system evolves through 
two distinct crossover temperatures into a very simple band structure with one free c-electron band 
and a disconnected upper and lower Hubbard band for /-electrons. The / electrons thus drop out 
of the Fermi surface volume as the temperature increases. 
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The electronic structure of heavy electron compounds 
remains a largely unresolved problem of solid state 
physics These materials show two characteristic 

'crossovers temperatures' in their transport properties, 
which arc commonly called the Kondo temperature, Tk , 
and the coherence temperature, T co h<TK- At temper- 
atures < T co h the materials are Fermi liquids, whereby 
the strongly correlated and almost localized /-electrons 
do participate in the Fermi surface volume. This has 
been established by a wide variety of de-Haas van-Alphen 
(dHvA) experiments That the Fermi surface vol- 

ume in heavy electron metals may be subject to dras- 
tic changes, however, has been demonstrated by dHvA 
experiments in the neighborhood of the so-called meta- 
magnetic transition Q. This transition appears to cor- 
respond to the /-electrons 'dropping out' of the Fermi 
surface volume. While the requirement of long quasi- 
particle lifetime restricts dHvA experiments to very low 
temperatures and thus does not usually allow to scan the 
temperature evolution of the Fermi surface, it has been 
conjectured |(| that the /-electrons do also 'drop out' of 
the Fermi surface volume as the temperature increases 
above T co h- The observation that the iso-structural com- 
pounds CeRu2Si2 and CeRu2Ge2 have Fermi surfaces, 
which differ rather precisely by the volume associated 
with the single /-electron/unit cell ||J|,[jJ would be quite 
consistent with this scenario if one were to assume |||7]] 
that T co h is higher than the experimental temperature 
for CeRu2Si2 but lower for CeRu2Ge2. More precise in- 
formation about the temperature-dependence of the elec- 
tronic structure of the Kondo lattice therefore would be 
desirable. In the following we want to present informa- 
tion from computer simulations. We have performed a 
Quantum Monte Carlo (QMC) study for a simple ID 
tight-binding version of the Kondo lattice (or periodic 
Anderson model), with Hamiltonian 

H = -*5>U«rV + H.c.) -Vj2(4,Ji„ + H.c) 



Here cj a (f] a ) creates a conduction electron (/-electron) 
in cell i, n i ^=f} (T f i a . In the case of 'half-filling' i.e. two 
electrons/unit cell and in the so-called symmetric case, 
ef=U/2, the model has particle-hole symmetry whence 
the minus-sign problem of the QMC procedure is absent 
and very low temperatures can be reached 0. Surpris- 
ingly enough we found that this highly desirable feature 
persists to a high degree also away from half filling. For 
example even at temperature T—0.05t, corresponding to 
1.2% of the conduction electron bandwidth, we obtained 
a QMC expectation value for the Fermionic sign of 0.605 
and thus obtain very well-converged spectra. This al- 
lowed to perform a detailed temperature scan of the en- 
tire dynamics of the model. In the following, we use 
the parameters U=2ef=8t, V=t. Throughout we have 
fixed the chemical potential at ji=—t\ this gives an elec- 
tron density of n—1.6, with variations of less than 1% as 
the temperature changes from 0.03i to t. The 'all elec- 
tron Fermi momentum' Uf then is 0.8-7T, the 'frozen core 
Fermi momentum' obtained for unhybridized conduction 
electrons (i.e. V=0) would be fc^=0.37r. 
Figure |l| shows the single particle spectral function at the 
lowest temperature for which it was calculated, T=0.05i. 
Thereby two different spectra are shown, one of them for 
electron creation/annihilation in the c-orbitals, the other 
one for /-orbitals. The spectra contain several distinct 
features, which are consistent with previous exact diag- 
onalization |)| and analytical [jl0| results: at the high- 
est excitation energies, we have two dispersionless /-like 
Hubbard bands, separated by an energy of rs U. In be- 
tween these, a c-like tight-binding band with dispersion 
£fe=— 2t cos(fc) can be seen. This band seems to cross the 
chemical potential close to k F . Closer inspection of the 
spectra shows, however, that this is not the true Fermi 
momentum of the system. Very close to the chemical 
potential (see Figure for a close-up of 
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FIG. 1. Single particle spectral function for the 16-unit cell 
system. 

this region) there is a rather intricate network of /-like 
bands with low spectral weight. Some of these bands 
apparently hybridize with the c-like band. It should be 
noted that the peak positions of the c- and /-like band 
do not always agree perfectly but actually the MaxEnt 
procedure reconstructs the spectral densities of the two 
species separately and does not use the peak positions 
of the other species as input information. The possi- 
ble discrepancies between c-like and /-like peak positions 
thus give a feeling for the accuracy of the results. For 
most momenta a total of three /-like bands can be dis- 
tinguished. In the outer half of the Brillouin zone there 
is one band right at /x, which seems to disperse through 
fi approximately at kp, the Fermi momentum for c and 
/-electrons combined. This band seems to hybridize with 
the c-like cosine-band, producing the familiar hybridiza- 
tion band structure indicated in the Figure. In addition 
there are two /-like sidebands at slightly higher excita- 
tion energies ~ ±0.5i. While we cannot rule out that 
these side-bands do have a tiny dispersion, they never 
approach the chemical potential. The three-band struc- 
ture around /x can also be seen very clearly in the angle- 
integrated /-electron spectral density, see Figure |[ The 
angle-integrated spectrum is actually not obtained by 
averaging over the angle resolved spectra. Rather this 
quantity is obtained by an independent measurement and 
MaxEnt run, and thus might be viewed as an indepen- 
dent cross check for the band structure. All in all the 
low energy band structure is quite reminiscent of our 
earlier results for the Kondo insulator particularly 
the /-like side-bands could also be clearly identified in 
this case. With the exception of the /-like sidebands one 
can describe the low-energy band structure phcnomcno- 
logically by the hybridization of the c-like cosine band 
with a practically dispersionless 'effective' /-level, which 
is pinned to the chemical potential of the unhybridized 
conduction electrons. This is precisely the picture in- 
ferred from previous exact diagonalization || and ana- 
lytical |l(| calculations. 

We proceed to a discussion of the temperature depen- 



dence of the band structure around /x. As the temper- 
ature is increased to T=0.1t, the /-like bands near /x 
generally become quite broad. Near fc=7r the upper 
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FIG. 2. Single particle spectral function at different tem- 
peratures. The thick dotted line for T=0.05t shows a plausible 
form of the low-temperature band structure. 

side-band becomes very diffuse and shifted toward higher 
energies 1.5t. The lower side band still can be clearly 
distinguished, and also near fc=0 an /-like part slightly 
above /x is clearly resolved. This can also be seen in 
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the angle-integrated /-density (Figure ||) : the lower side- 
band and the band at /i persist, but the upper side-band 
is shifted to relative large energies. In between these two 
bands there is still an indication of the dispersive /-band 
which crosses fi. This changes completely at T—0.2t, 
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FIG. 3. Angle-integrated /-electron spectral weight. 

where the band crossing through // has disappeared and 
only the side-bands remains. The lower one seems to 
have shifted to somewhat lower excitation energy, but 
does not reach /i, see also Figure ||. The only dispersive 
band which could intersect [L now is the c-like conduction 
band, which crosses n near kp. The low-energy /-like 
band which intersected n at /cf has disappeared, and 
the only /-like feature near /i are the two side-bands. 
These are practically dispersionless and never even ap- 
proach /i, so that there is no more band which could cross 
u at the 'all-electron Luttinger momentum' kp. This 
suggests that at this temperature the Fermi surface is 
purely c-like, with the smaller Fermi momentum k F : the 
/-electrons have dropped out of the Fermi surface vol- 
ume. Precisely the same behavior has previously been 
observed in the Kondo insulator: upon increasing T the 
/-like bands which formed the edge of the insulating gap 
at low temperatures disappeared, the dispersionless side- 
bands persist, and the Fermi surface becomes c-like. The 
only difference is that in the case of the Kondo insulator 
the collapse of the Fermi surface volume is equivalent to a 
insulator-to-metal transition Q . Proceeding to the high- 
est temperature we have studied, T=t, we note that there 
even the two /-like sidebands have disappeared and the 
Fermi surface now is obviously purely c-like, with Fermi 
momentum k F . The angle-integrated /-electron density 
shows no more spectral weight at fi, the only remain- 
ing /-like feature on the photoemission side is the lower 
Hubbard band at energy w — 3.2i. Increasing the tem- 
perature thus leads to a transfer of /-like weight from the 
Fermi energy to the Hubbard band. With the exception 
of the side bands the overall shape of the /-density at 
low temperatures is also quite consistent with the results 
from 1/N expansion for the Anderson impurity model 
p2[ . Also, the overall evolution with temperature, in- 
cluding the disappearance of all /-like weight from the 
region around /i at high temperature is completely con- 
sistent with the situation for the Kondo insulator [§| . 



To get a more detailed picture of the T-dependence of 
the Fermiology, we consider the momentum distribution 
n Q (fc)=^ .(a^. cr afc cr ), where a—c, /. We focus on the /- 
like distribution and start out at the lowest temperature. 
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FIG. 4. Momentum distributions at different tempera- 
tures. The thin dashed lines give the Fermi-Dirac distri- 
bution computed with eu = — 2£cos(fc) and y, = — t at 
T = t,0.5t,0.2t. 

As k increases from the center of the Brillouin zone, 
there is first of all a steep rise of n/ (k) at approximately 
k F . Inspection of the single particle spectrum (Figure 
||) reveals that this is due to an increase of the spec- 
tral weight in the side-band for to — fi < - note that 
nf(k)= J du;Af(k,LLi)f(uj), where f{ui) is the Fermi func- 
tion. As we approach k=n and pass through kF the 
/-distribution drops again, which this time is associated 
with the crossing of the central /-like band through fi 
(see Figure |||) . Due to the flatness of this central band 
and its small spectral weight per fc-point, Z, the 'drop' is 
actually very shallow; in a Fermi liquid the momentum 
distribution near kp can be written as 

. . Zvp . . 

n[k) = const + (k — kp), 

k B T 

where Z is the quasiparticle weight, vp the Fermi veloc- 
ity and const the contribution from the incoherent back- 
ground (assumed to be ^-independent). Both, a small Z 
and a low Fermi velocity will therefore make the slope 
of n{k) very small. However, for the lowest tempera- 
ture, T = 0.033£, one can envisage how the Fermi step 
sharpens up at low temperatures. At T=0.1t, the drop 
around kp is still visible, but becomes very weak. Such a 
flattening of the Fermi edge discontinuity with increasing 
temperature is quite natural. Inasmuch as the discrete 
/c-mesh allows for such a conclusion, the four nf(k) for 
T < 0.1 all intersect approximately at kF - this is what 
one would expect if kF were the Fermi momentum (on 
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the other hand, the same holds true also for k F , where 
inspection of Af(k,uj) suggests that the drop is associ- 
ated with a reduced spectral weight of the side band). 
At the intermediate temperature T=0.2t no more drop 
of rif(k) can be distinguished near /cf, and at T=t the /- 
distribution becomes completely flat indicating the more 
or less complete localization of the /-electrons. 
Turning to the c-electron distribution we first of all note 
that all curves intersect at approximately k F , and that 
the common value of nf(k) is unity at this momentum. 
This would suggest that k F is the Fermi momentum in 
a Fermi liquid with Z—l and for T > 0.5i n c (k) is in 
fact nearly identical to that of free, unhybridized con- 
duction electrons. This changes as the temperature falls 
to T=Q.2t, where n c {k) becomes temperature indepen- 
dent. This would be consistent with the assumption that 
the formation of the heavy band below T=0.2i the true 
Fermi surface is shifted from k F to kp, so that below 
this temperature the drop of n c (k) at k F is due to k- 
dependent hybridization with the heavy /-band, rather 
than a Fermi surface crossing. 

Taken together, the data for the single particle spectral 
function and the momentum distribution rather clearly 
suggest a temperature induced volume collapse of the 
Fermi surface. At low temperatures a nearly disper- 
sionless band of low /-like spectral weight crosses the 
chemical potential near /cp, the Fermi momentum for c- 
electrons and /-electrons combined. Thereby in agree- 
ment with analytical work JlCfl the heavy band is pinned 
to the chemical potential of the Fermi sea of unhybridized 
conduction electrons. The momentum distribution of the 
c-electrons still shows a sharp drop at k F [ pOlpl] ] which 
becomes T-independent at low T and thus cannot develop 
into a Fermi edge at T—0. As the temperature increases 
the heavy band disappears, whence the Fermi momen- 
tum must collapse to k F . Consequently the /-electron 
momentum distribution then shows no more notable vari- 
ation near kp, whereas the drop of the c-electron distri- 
bution near k F becomes temperature dependent. 
In summary, we have studied the single particle spectrum 
of the one-dimensional Kondo lattice. A spectacular fea- 
ture from the technical point of view is that the notorious 
minus-sign problem is nearly absent in this system, so 
that reliable QMC simulations can be performed even at 
very low temperatures. This opens the possibility to per- 
form detailed numerical studies of the nontrivial temper- 
ature evolution of its electronic structure. Thereby full 
advantage can be made of the special features of QMC, 
such as the momentum-resolved calculation of all dynam- 
ical correlation functions. 

Our data show that despite its one-dimensional char- 
acter, which would place the fully interacting system 
into the universality class of Luttinger liquids, the sys- 
tem shows a behavior which is quite reminiscent of the 
Fermi liquid-like three-dimensional Heavy Fermion com- 



pounds. In particular we could find the equivalents of the 
two characteristic temperatures of Heavy Fermions, the 
Kondo and the coherence temperature. Both characteris- 
tic temperatures are associated with a distinct change of 
the electronic structure near /i: at the lower temperature, 
which we associate with the coherence temperature, the 
'heavy' /-like band, which crosses the Fermi energy and 
forms the Fermi surface, disappears. Above this temper- 
ature the Fermi surface is c-like, and does not comprise 
the /-electrons any more. Only two dispersionless /-like 
sidebands remain at relatively low but finite excitation 
energy. Around the second characteristic temperature 
these side bands disappear as well, so that the c- and 
/-electron system become more or less decoupled. 
In our previous work on the Kondo insulator we have 
also shown that the temperature where the /-electrons 
'drop out' from the Fermi surface volume depends sensi- 
tively on the c-f hybridization V (this is also to be ex- 
pected on the basis of the impurity model fll^] ), and we 
may expect that the same holds true also for the metal- 
lic case. At a given temperature different compounds 
with different c-f hybridization, such as CeRu2Si2 and 
CeRu2Ge2, may therefore be in different regimes, whence 
their Fermi surface volume would differ by the number 
of /-electrons. 
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